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Abstract 

This paper deals with a single-server queue with modulated arrivals, service requirements 
and service capacity. In our first result, we derive the mean of the total workload assum¬ 
ing generally distributed service requirements and any service discipline which does not 
depend on the modulating environment. We then show that the workload is exponentially 
distributed under heavy-traffic scaling. In our second result, we focus on the discrimina¬ 
tory processor sharing (DPS) discipline. Assuming exponential, class-dependent service 
requirements, we show that the joint distribution of the queue lengths of different cus¬ 
tomer classes under DPS undergoes a state-space collapse when subject to heavy-traffic 
scaling. That is, the limiting distribution of the queue length vector is shown to be expo¬ 
nential, times a deterministic vector. The distribution of the scaled workload, as derived 
for general service disciplines, is a key quantity in the proof of the state-space collapse. 
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1 Introduction 

Markov-modulation is a way to formalize the embedding of queues in a random environment. 
The parameters of the queue in question, typically arrival rates, service requirements or both, 
are governed by an external Markov chain, thereby creating an extra layer of randomness 
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around the stochastic queueing process. For classical resnlts on Markov-modulated single¬ 
server queues with the first-come-first-serve (FCFS) discipline see e.g. [271 HOj [3T] . Recent 
work on systems in a Markov-modulated environment can be found in for example dZlIMlES]. 
In this paper we will analyse a modulated queue under a heavy-traffic scaling, that is, evaluate 
the system at its critical load. It is a well-known result from the literature of single-server 
quenes, that under fairly general conditions |23j, the steady-state distributions of the ap¬ 
propriately scaled queue length and workload become exponential when the critical load is 
approached. This property has been seen to carry over to certain systems where arrivals 
and service times are modulated by an external Markov process, see o HSl dS]. In fact, [3] 
establishes an even stronger result: convergence of the queue length process to a reflected 
Brownian motion. Multi-class single-server queues under a heavy-traffic scaling have been 
studied in e.g. [22j for FCFS with feedback routing, [20] for the discriminatory random order 
of service discipline, and [141 [33] for the discriminatory processor sharing (DPS) policy. In 
particular, mm show that the steady-state queue length vector undergoes a so-called state- 
space collapse and converges to an exponentially distributed variable, times a deterministic 
vector. The cited multi-class results under heavy-traffic scaling are all for non-modulated 
systems. In light of this, we will in this paper put special emphasis on a modulated multi-class 
single-server queue, and the limiting steady-state queue length distribution is derived. 

While there is little ambiguity in how arrival rates are modulated, there are in the literature 
typically two ways in which to modulate the service rates. One can (i) let the departure rate 
be continuously modulated throughout a customer’s service, the other approach is to (ii) let 
a customer’s service requirement distribution be based on the state of the environment when 
it arrives and remain the same until the customer departs. We note that by adapting the 
number of different customer classes, the fixed service requirements of case (ii) can be seen 
as a special case of the continuously modulated requirements (i); we further elaborate on this 
later in the paper in Remark [1] in Section [2l 

The way the load or traffic intensity for modulated queues is defined goes hand in hand 
with the way the service rates are affected by the environment. In case (i), the load is 
typically the average of arrival rates (where the averaging is with respect to the equilibrium 
distribution of the environment), say Aqo, divided by the average of service rates, say Hoo 
(see e.g. [27] 1. In case (ii), the load is taken as the average over the arrival rate times the 
mean service requirement, say, Xd/g-d per state d of the environment (see e.g. [Ill [31]). The 
two load definitions represent different scenarios. In particular, when load (ii) is equal to 1 
(the critical load) it means that in at least one state of the environment, the total load over 
all classes must exceed 1, i.e. for at least one state we must have overload. This is true only 
for special cases of definition (i). 

In this paper, special focus will be given to a multi-class single-server queue under the DPS 
discipline. The DPS discipline was first introduced by Kleinrock in [25] as an extension of 
the well-known egalitarian processor sharing (PS) discipline and has turned out to be very 
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suitable to model the simultaneous parsing of diverse tasks, such as processing network data. 
Under this service discipline, the service capacity is divided between all present customers in 
proportion to their prescribed weights. Due to the challenging nature of DPS systems, most 
available results are in terms of limit theorems and moments. Fayolle et al. m established 
the mean sojourn time conditioned on the service requirement. That analysis also yielded 
the mean queue lengths of the different classes, which were shown to depend on the entire 
service requirement distributions of all classes. The DPS model has finite mean queue lengths 
irrespective of any higher-order characteristics of the service distribution, see Avrachenkov 
et al. [S]. This is an extension of a result for the Processor Sharing (PS) system, which holds 
while the queue is stable. DPS under a heavy-traffic regime was analysed in [TJ] assuming 
finite second moments of the service requirement distributions. Assuming exponential service 
requirement distributions, a direct approach to establish a heavy-traffic limit for the joint 
queue length distribution was described in m and extended to phase-type distributions 
in [33] . Combining light and heavy-traffic limits, in m an interpolation approximation is 
derived for the steady-state distribution of the queue length and waiting time of DPS. The 
performance of DPS in overload is considered in [2]. Asymptotics of the sojourn time have 
received attention in UM- Game-theoretic aspects of DPS have been studied in [3H [18]. A 
thorough overview of DPS results can be found in [T]. 

We are not aware of work analysing a DPS system under modulation. We refer to [28] 
where the Processor Sharing discipline (DPS discipline with unit weights) was analysed in 
a Markovian random environment. Multi-class queues in a random environment have been 
studied for different models in [II1[32|. In [32] , a modulated system is studied where arrivals 
can only occur at transition epochs of the modulating process but service requirements are 
class-dependent and generally distributed. Using a time-changing argument, the waiting 
time distribution is derived under the FCFS discipline. In the authors derive a Brow¬ 
nian control problem to establish a form of the c/x scheduling rule in heavy traffic under 
continuously modulated service requirements. By using a particular scaling, the time-scale 
separation of the external environment and the queue length process is exploited. Similar 
scaling of a modulated queue can also be seen in results on the Markov-modulated infinite 
server queue in e.g. [7]. 

The system we analyse in this paper is a single-server queue where the arrival rates, service 
requirements and service capacity are modulated. We focus on the setting where a customer’s 
service requirement distribution is based on the state of the environment when it arrives 
and does not change throughout its service. The service capacity is however continuously 
modulated. This assumption is in line with the literature for various types of modulated 
queues, see [iniiiai2siE2]. In Remark |2| in Section [U] we discuss how part of our results can 
be extended to a more general model with continuously modulated service requirements. We 
derive the distribution of the workload under a heavy-traffic scaling for generally distributed 
service requirements and any service discipline which does not depend on the environment. 
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We then turn our attention to the DPS discipline in a multi-class queue, which is a particular 
case of the general modulated system as described above. The weights of the DPS system 
determining the service proportion, depend on a customer’s class and do not change with 
the environment, which means that the workload result remains valid. An important finding 
in the present paper is that the queue length vector under DPS becomes independent of the 
modulating process in the heavy-traffic limit, which is consistent with the modulated M/G/1 
queue in e.g. mm- This, together with the obtained result on the workload, allows us to 
derive the distribution of the queue length vector under DPS and to show that it undergoes 
a state-space collapse. 

The remainder of the paper is organized as follows. In Section [2] we describe the model. In 
Section [3] we study the workload of a single-server queue with modulated arrivals, service 
capacity and service requirement distribution, establishing its distribution in heavy traffic. 
From Section m onwards, the focus is on DPS. In Section 0] we derive some basic properties 
of the queue length distribution, obtain a rate conservation law and derive an equation for 
the moments of the queue lengths weighted with the modulated service capacity. Section [5] is 
devoted to the heavy-traffic scaling; there we show that the distribution of the environment 
becomes independent of that of the queue length vector, in addition to deriving two technical 
lemmas. The exponential limiting distribution of the joint queue length in heavy traffic 
follows in Section [6l The result is shown in two steps in the subsections 16.11 about the state- 
space collapse and 16.21 about the exact limiting distribution, where we rely on the workload 
result of Section [3l We conclude with a summary and some open questions in Section [71 

2 Model 

We analyse a single-server queue modulated by an independent external environment, which 
is formalized by an irreducible continuous time Markov chain on a finite state-space {1,..., D}. 
The modulating process is denoted with Z and is governed by an infinitesimal generator ma¬ 
trix Q = {qd£)di=i with an invariant distribution tt = (tti, ... In what follows, vectors 

are generally denoted in bold. New customers arrive according to a Poisson distribution 
with rate Xd when the environment is in state d. A customer arriving in state d has a service 
requirement distribution given by a function Hd{-) with Laplace-Stieltjes transform (LST) 
hd{-)- The service requirement does not change further with the environment. The first and 
second moment are given by hdi and hd 2 , respectively. In addition we let the service capacity 
be scaled by a factor Cd during the environment’s stay in state d, this can thus change during 
the service of the customers. The traffic intensity will be measured as 

Poo — ^ ^^ '^dXdhdl/C OO; (1) 

d 

with Coo := T^dCd being the service capacity averaged over the environment. The workload 
is defined as the time it takes to empty the system at an arbitrary moment in time given 
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the observed environment and is denoted hy W. In Section [3] we study the workload and 
the environment as a two-dimensional process {W,Z), under any service discipline that is 
independent of the environment. 

The first main result of this paper concerns the distribution of the workload when the traffic 
intensity approaches its critical point. The system is said to be in heavy traffic when p^o 
approaches 1. Let > 0 and define the following parametrization 

^ ^d, as iV-)■ oo, (2) 

poo poo 

where poo is based on the unsealed parameters. Prelimit quantities will be denoted with a 
superscript (N)] the prelimit traffic intensity is thus and is equal to 1 — 1/N. Limiting 
quantities will have a in heavy traffic the traffic intensity is denoted poo and is equal to 1. 
In the remainder of the paper, starting in Section 01 we analyse a single-server queue with 
K customer classes under the discriminatory processor sharing policy; the queue is again 
embedded in a random environment. Let be the probability that a customer, that 
arrives while the environment is in state d, is of class k; note that a^^d = 1 for a given 
d. The Poisson arrival rate of a class-A: customer is denoted with X^^d ■= (^k,d^d and for each 
class k it is assumed that Xk^ > 0 for at least one state d. In the multi-class setting we assume 
that a class-A: customer has an exponentially distributed service requirement with mean 1/pk- 
We believe that the results obtained in this paper can be extended to phase-type distributed 
service requirements, the latter being dense in the space of all distributions on [0, oo). For 
the non-modulated DPS queue, the phase-type analysis was performed in [33] using similar 
proof techniques. For ease of exposition, however, we focus here on the exponential case. 

We no longer let the service requirement of a particular customer be environment-dependent 
(although the distribution of an arbitrary customer is, as explained in Section 16.2p . One can 
however retrieve the environment-dependent service requirements by introducing additional 
classes for each environment, see Remark [TJ By referring to a class-A; customer’s service 
rate while in state d as pk,d '■= d’k^d-, we take the modulated service capacity into account. 
Most of the results for the queue length can in fact be shown without assuming this product 
form, representing a system where the service requirements are continuously modulated, see 
Remark [31 Section f6.ll for further details. 

We denote the average arrival rate of class-A; customers by Xk^oo '■= '^d^k,d'^d, similarly 
we denote the average service rate for class k by Pk,oo ■= '^dk-k,d'^d = d-kC-oo and Pk,oo '■= 
^k,oo/d‘k,oo- The aggregate traffic intensity for the multi-class model is defined as 

K 

Poo ■ ^ ^ Pk,oo) 

k=l 

which is consistent with the definition in Eqn. ©• 

Let the state of the multi-class system be described by the vector of random variables 
(Ml,..., Mk, Z) =: (iW, Z), where Mk is the number of class-A: customers, for A: = 1,... , iL. 
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As before, Z represents the state of the background process. In a DPS system, the random 
fraction of service given to a class-A: customer is 

9k 

9j^Ij 

where gk are weight parameters associated with each class k. 


When in heavy traffic, we denote pk,oo '■= ^k,oo/9k,oo and thus have for the multi-class system 


,CXD 1 Afc ,oo 1 

Pk, QO — / ^ — / , — / , Pk, QO — I- 

^ ^ 9k,oo Poo^Pk,oo Poo^ 


Remark 1 (Modulated service requirements rewritten to classes). Any multi-class system 
where the service requirement distribution is determined by the modulating process at a 
customer’s arrival, can be written as a multi-class model with non-modulated, only class- 
dependent, service rates as illustrated below: While in state d of the environment, class-A: 
customers arrive with rate and have exponential service requirement with mean l/pk,d 
and weight pk- Such customers we refer to as being of class {k,d) and count with M^^d, 
hence we need to keep track oi K ■ D different customer “classes”. Arrivals to class (A:, d) are 
inactive while not in state d. 


classes 

arrival rates 

serv. rate 

weight 


d = l 

d = 2 

de {1,2} 

(1,1) 

Ai,i 

0 

9i,i 

91 

(1,2) 

0 

Ai,2 

9l,2 

91 

(2,1) 

^2,1 

0 

92,1 

92 

(2,2) 

0 

^2,2 

92,2 

92 


Table 1: A multi-class system where the service requirement distribution is fixed upon arrival 
can be translated into one with non-modulated service requirements. 


From Table [T] one can easily see how a K = D = 2 system can be written into one with 
A' = 4 and D = 2. The arrival rates are still modulated, but in an on-off way. The service 
rates pk,d are now non-modulated. 


3 Workload 

In this section we consider the workload in a modulated queue and extend the results for an 
M/G/1 type queue from Falin and Falin [15] and Dimitrov [T2| to include modulated service 
capacity. We derive the mean of the workload and then its distribution in the heavy-traffic 
regime. 
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Let po,d = P{W = 0, Z = d) and let a = (oi,... be a vector solving 

[Q ■ ^]d — Cci ^dhdl Coo(l Poo)) (3) 

for d = 1,... ,D. Note that such a solution always exists since the right hand side vector of 
Eqn. ([3]) is orthogonal to tt. We obtain the following result for the mean workload. 

Proposition 3.1. For any service requirement distribution and any service discipline 

that does not depend on the state of the environment, the mean of the workload satisfies 

Coo(l Poo) 

where a is a solution of Eqn. ([3]). Furthermore, 1 — poo = X^dPo.rfCd/coo- 

Remark 2. Although the solution vector a is not unique, the term 

D 

X y C(i) A Po,d*"d]) 

d=l 

as appearing in Eqn. (j3|), is. This is due to the following argument: Suppose a and a* are 
two solutions to Eqn. ([3]). Then 0 = Qa — Qa* = Q{a — a*), so {a — a*) is in the nullspace 
of Q. But Q is a generator matrix so it can easily be seen that Qr = 0 for any vector 
r • 1^ = (r, r,... , r)^, r G M. Also, since the environment is an irreducible Markov chain, the 
nullspace of Q has dimension 1, and therefore, (a — a*) = ra ■ 1^, for some G M. Thus, 

D 

'^{ad - a*^)[7rd{Xdhdl - Cd) + P0,dCd] = raCoo{Poo - 1) + raCoo(l - Poo) = 0, 

d=l 

where the hrst term follows by definition of poo and Coo and the second term comes from 
Eqn. ([7]) in the proof below. 

Proof of Provosition \d.l\. Define Fd{x,t) = P{W{t) < x,Z{t) = d), for some time t > 0. 
In an infinitesimal time df, a new arrival requiring service x changes the workload with 
probability XdPfd{x)dt. The service capacity is scaled by Cd when the environment is in 
state d, meaning that in dt time, the workload is reduced by c^dt, yielding by a classic 
birth-and-death argument for the M/G/1 queue. 


Fdix, t + dt) = {l- Xddt + qdddt)Fd{x + Cddt, t) 

pX 

+ y] qidPdx + Cidt, t)dt + Xddt / Fdix + Cddt - y, t)dHdiy). 

Jo 


e^d 


We let t —)• oo to go to steady-state and since 


Fdix + Cddt) - Fdix) Fdix + Cddt) - Fdix) , 

-;-= Cd -;- —^ Cdkj(x) 

dt Cddt df4,o 
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we obtain 


CdF'di^) = i^d - qdd)Fd{x) - ^ qidFiix) -Xd Fd{x - y)dHd{y). 

ij^d 

Denote the LST of Fd{-) by ipdis) = Z{t) = d] = po^d + /o+6~^*dFrf(x). The 

corresponding transform equation becomes 
D 

^ qidq^iis) = [Arf(l - hd{s)) - scd](pd{s) + spo,dCd. (5) 

i=i 

It is now convenient to sum over d and divide through Eqn. ([5]) with s to get zero on the 
left hand side, leading to 

(1 - hd{s)) 


d d 

Using that (pd{0) = T^d and by I’Hopital 

(1 - hd{s)) 


Cd - Xd- 


lim '' = — lim h'j(s) = hd ^, 

sto s sto 


we get by taking the limit s ^ 0 of Eqn. ([6|) that 

P0,dCd 


E 


- 1 Poo- 


( 6 ) 


(7) 


We differentiate Eqn. ([6|) w.r.t. s: 

J2^dis)Xd 


h'di^) ^-hd{s) 

' O 




Cd - Xd 


1 - hd{s) 


which in the limit of s —>■ 0 results in 




Coo(l poo) ^ ^ d — ^ ^ 


H“ ^d'^di.^dhdl Cc/) “1“ P0,dCd^d 


d d 

which gives the desired expression for the mean of the workload, Wd- 


( 8 ) 


d d 

with the first moment of the workload while in state d being 

- lim = E[W^, Z = d]=: Wd. 

siO 

Now multiply Eqn. ([5|) with ad, sum over d, take the derivative w.r.t. s and let s —>■ 0 to 
obtain 

^ ^ Wd[Coo(Poo 1) T Crf Xdhdl] — ^ ^ [cdT^diXdhdl Crf) T Po,dCdCd\ ; 
d d 

by using Eqn. Adding this equation to Eqn. ([5]) yields 

'^dXdhd2 


( 9 ) 

□ 


















Eqn. ([7|) makes it clear that in heavy traffic, that is when —>■ 1, all the probabilities 

p^J go to zero. Also, in heavy traffic, the right hand side of Eqn. ([3]) reduces to — Xdhdi- 
Recalling the parametrization 1 — = 1/N in Section [21 we obtain the following result. 


Proposition 3.2. For any service distribution Hd{-) and any service discipline that does not 
depend on the state of the environment, the mean of the workload in heavy traffic satisfies 


lim 

TV— >-CX) N 



^dhd2l^ ^di^dhdl ^d) ; 


( 10 ) 


where a is a solution of [Q ■ a\d = Cd — Xdhdi- 

Proof. Under the heavy traffic scaling, the empty probabilities p^J go to zero for d = 
1,... ,D. The result then follows immediately from Eqn. Q and 1 — = 1/N. □ 


This leads to the main result of this section. 


Theorem 3.3. In heavy traffic, the scaled workload , converges in distribution to W, 

where W is exponentially distributed with mean given in Eqn. m- 

Proof. This follows from combining Proposition 13.21 with Theorem 4 in [12]. The full proof 
is in the Appendix. □ 

The above results yields that W is relatively compact, which together with the metric space 
being separable and complete implies that the scaled workload is tight, by Prohorov’s 

theorem [6|. 


4 Queue length vector under DPS 


In the remainder of the paper we focus on the multi-class model under DPS, where we assume 
that the service requirements of class-A: customers are exponentially distributed with rate pk- 
In this section we establish some properties of the joint queue length distribution. We start 
with the flow equations, followed by a rate conservation law and an equation for the moments 
of the queue lengths conditioned on the environment. 

Equating the flow in and out of state {M,Z) = {m,d) yields (noting that —qdd = 


f 'y ^X Xk,d + Pk,d ' l{rrn.>0}) Idd j Pm,d 

\k=l / 

y yXf^ dPm.—ep;,d ■ l{mu>0} ~k i Pk,dPm-\-ei.,d) T ^ ' QdlPm,lj (H) 

^ Hihimi + gk 
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where Pm,d ■= = {m,d)) and is the vector with 1 in the k-th. place and zeros 

elsewhere. We now define the partial probability generating function (PGF) for when the 
background process is in state d: 


Pd{z) ■hz=d}\ 

oo oo 

^ ^ ¥{Mi = mi,... ,Mk = mK,Z = d) ■ 

mi=0 

= Y, Prr.,dZ^, 

m>0 


where z^^ ■ ■ ■ z^^ =: and (mi,..., mx) > (0; • • • j 0); be- m > 0. Then the overall 

generating function for the queue length is P{z) := • • • z^^] = X]£=i Pdi^)- We also 

dehne 


Rd{z) 


E 

m>0 


Pm,dZ'^ 

EjPjrrij 




hence 


dRdjz) 

dzk 


= z, 


E 


rrik 


m>ek 


Pm,dZ 


{EEi"ii>0}' 


By multiplying the flow equation (jlip with z'^, summing over all vectors m > 0 and rear¬ 
ranging terms, we can eventually write it in terms of the PGF Pd{z) and the partial derivative 
dRd{z)/dzk, that is. 


K 

E 

k=l 


Z}^P(i^z'j -\- fdk,d9k(.Zl^ 1 ) 


dRdjz) 

dzk 


D 


T.Pik‘ 


( 12 ) 


1=1 


It will be convenient to write the equation fully in terms of dRd/dzk, so we note the relation 

^ dRdjz) 
ykZk 

k=l 


D f \ dRdjz) 

Rdjz) = - \-po,d 


(13) 


where po,d = F((iW,Z) = j0,d)) is the probability of an empty queue in state d, which is 
equivalent to the probability of no workload dehned in Section [3l We incorporate this into 
Eqn. (fT^ to obtain 


K 


^ ^ ^fc,d(f Zff^) 


k=l 


K 


dRdjz) 

- +P0,d 


i=i 


dzi 


K fix f 1 ^ ^ 

+ Y P-k,d9k jzk - 1) ^ =Y '^dkZk 

k=l ^ 1=1 U=1 


dRjjz) 

dzk 


-Po,e 


Qed- 


(14) 


This equation we will use later when deriving the heavy-traffic limit. 

For the M/M/1 queue with modulated arrivals and service times, moments of the queue 
length and the sojourn times can be found for the FCFS service discipline, in [35] and |26j . 
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respectively. In [30] the authors establish a recursive formula to calculate moments of the 
queue length in a non-modulated DPS system. In a similar fashion, we obtain an expression 
for the sum of the state-dependent moments weighted with the capacity of the server. We 
also derive a rate conservation law, which shows how the average arrival rates per class are 
proportional to the resources allocated to that same class, and the service they receive. Both 
results can be found in the following proposition. 


Proposition 4.1. When the queue is stable, the average number of class-k arrivals is pro¬ 
portional to the service resources allocated to class-k customers, i.e. 


^k,oo ^ ^ hkjdM' 


Yhj 9jMj 


■ ■ l{Z=d} 


for k = I,..., K. Furthermore, the state-dependent expectations of Mk satisfy 

tuTA/T 1 I ^k,cx> , ^k,d^[Mj ■ l[z=d}] + ^j,dJ^[Mk ■ '^{Z=d}] 

2^ CdE[Mk ■ l{z=d}] = - + 2^ 9j - 


Fk 


Idk9k + k'j9j 


d d,j 

Proof. We sum Eqn. (|I2I) over d and take the derivative w.r.t. Zi to obtain 


E 


-K,dPd{z) + ^j,d{^ 


- Zi 


dPdiz) dRd{z) 


Letting z —>■ I yields 


Since 


E 


d 




dzi 


dRd{z) 




dzi 


'^Fj,d9j{zj - 1) 


d’^Rdjz) 

dzidzj 


dzi 


^i,dPd(,z) 


z^l 


z^l. 


= 0 . 


lim Pd{z) = TTd, 

z—^1 

the following conservation law results from Eqn. (|I5h . for k = 1,..., K: 

dRd{z)_ 

d “ " 


'^k,oo — / ^ fdk^d9k~ 


hk,d)^ 


dzk 

9k^^k 

22 j 9jMj 


■ l{E, M,>o} • l{z=d} 


= 0 . 


(15) 


(16) 


By taking partial derivatives of Eqn. (USD, we obtain after standard calculations the recursive 
relation 


d^Pd{z) 


dzi^ ■ ■ ■ dzi. 


K 


Y.9k 


d^^^Rd2 


z->-l k=l 


dzi^ ■ ■ ■ dzi.dzk 


^ • • • dzi 


(17) 


£=1 


In particular, this yields the explicit form 

dPd 


lE[Mfc • l{z=d}] = 


dzk 


z->-l 




d^Rd 


dzkdzj 


+ 9k 


dRd 


z^l 


dzk 


z^l 
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Proceeding from the rate conservation law, Eqn. dUD, and by using = ^kCd, we have 


E 


9-Rj 


^dQk 


dzk 




z^l 


l^k 


By taking two partial derivatives of the balance equation Eqn. (USD we can solve for a second 
mixed derivative of R^, namely 

d'^Rd _ ■ l{z=e}\qed + ^k,d)^[Mj ■ l{z=d}] + ^j,d)^[Mk ■ l{z=d}] 

z^l ^^k,d9k + 9j,d9j 


dzkdzj 


thus also yielding a mixed moment. Summing over the weighted moments of the number of 
class-A: customers while in state d, we obtain a linear equation resembling Eqn. (16) in [SOj 
for the non-modulated DPS queue: 


y^CrfIE[Mfc • l{z=d}\ 


^k,c 


l^k 


+ Y1 ^d9j 


d^Rd 


d,j 


dzkdzj 


(by Eqn. (fTTD ) 


z^l 


Afc,oo ■ l{z=e}\qid + Afc,dE[Mj • l{z=d}] + Aj,<iE[Mfc • l{z=d}] 

+ 9j - 


l^k 

^k,c< 




+ ^9j 


d,j 


Idk9k 9j9j 

Afc,dE[Mj • l|2=rf}] + Aj^rfE[Mfc • l^z=d}\ 

9k9k 9-j9j 


where the last equality comes from '^^qed = 0. 


□ 


5 Preliminary results for the queue length in heavy traffic 


We proceed to show that in heavy traffic, the distribution of the environment and the joint 
queue length become independent. This result, along with two technical lemmas that we de¬ 
rive in this section, will later help to establish the main result about the limiting queue length 
under DPS, presented in Section [6l Here we consider the queue length vector (Mi,..., Mk) 
scaled with 1/N and evaluate the PGP in . The objective is to determine the distribu¬ 
tion of ^(Mp\ ..., M^'^) ■ l^z=d} s-s goes to inhnity. We will state the existence of the 
limiting vector, and thus also the limit of the generating function ^^^( 2 ;^/'^), as an assump¬ 
tion. This assumption will be proven in Section [6.21 The superscript N denotes dependency 
on the prelimit parameter A^^\ 

We make use of the change of variables e~^'^ = Zk, for Sk > 0, and denote e“®/^ := 
(e-®i/'^,..., Assuming that the limit exists, we use the new variables to dehne 

the heavy-traffic quantities: We let po^d ■= ^ ^di^) ■= hmAr_).oo = 

limAr_^ooE[e“^j^^-’/^ • l{z=d}] and P{s) := Y^d^dis) = limAT-^oo E[e“^J We de¬ 

note by (Ml,... ,Mk) the random vector corresponding to the LST P{s). Finally, let 


Rd{s) := E 


Y,9jMj 


■ ^{E, Mi>o} • ^{z=d} 
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where the 1 in the numerator is to ensure that the bracketed expression remains bounded 
when the queue length quantities are all near zero. We can now proceed to the following 
lemma. 

Lemma 5.1. If liuiN^oo exists, then it satisfies 


K 


-Pd(s) = 


k=l 


dRdjs) 

dsk 


+ Po,d- 


Proof. From Eqn. (fT3l) . 


lim p{^\e = lim 


K 


N^oo 


N^oo ‘ 


9 k 


dRf\z) 


k=l 


dzh 


t=e-^/N 


+ PO,d- 


Note that 


dR^^'^(z) 

lim Zk —I- 

N->-oo OZk 


lim y 

m>e 

lim y 

W—>oo 


z=e~‘^/^ 

ruk 


m>e, 

ruk 


'Pm,d 


■ 1 




^—q—s/N 




(TV) ^-simi/N . Q—SKrTiK/N , ^ 


m 


m>ei^ ^ 


:Pm,d' 


lEj mj>0} 


= lim E 


N^oo 


M, 


(TV) 




(TV) 


e-E,-TMfViv . 1 




= E 


Mfc 




p .1 .. • 1 jT 

{T.jMj>0} ^{Z=d} 


dRdjs) 

dsk 


The second-to-last step follows from the fact that 


M 


-1 


(18) 


(19) 


( 20 ) 


{JliMjPyo} 


IS 


Ej 9j^j '-^J "'"j 

upper bounded by 1/ minjjgj). By the continuous mapping theorem (see Billingsley’s [6]), it 


converges in distribution to ^ • e ■ 1 rv- »> • The environment is not affected 

by the heavy-traffic scaling. Eqns. (jlOp and (j2np now conclude the proof. □ 

With the help of Lemma 15.11 we obtain: 

Proposition 5.2. //limjv^co T’j'^^(e“®/'^) exists, the joint queue length distribution is in¬ 
dependent of the environment in heavy traffic, that is 


Pdjs) = TrdPjs). 
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Proof. We use the change of variables = e • Since 


dRd, X 
Zk^—{z) 
dzk 


OSk 


we obtain, by applying the heavy traffic scaling to Eqn. oa, 




K 




dSn 

j = l J 

D 


{N) 

d 


f) 


E 


(e-/«)+pW 


E»- 


=1 Lfc=l 




g£d- 


With Taylor expansion we obtain 


E 


X (TV) / gfe Sfc 

' W iV2 


iV 


J=1 


D 

E 




K o„(TV) 


Sk , Sj^ \ oK^ 


fJ'k^dk ( jY 


9si 


( 21 ) 


=1 U=i 


(TV) 

.1 


qid + 0(iV 


-3^ 


Since —§f .— is bounded (see proof of Lemma 1.5,11) and converges to we obtain as 


(iV) 


N —>■ oo, 


D 


V ■ [Q]d = Y1 


e=i ik=i 


K 


dRAs) 

2^gk ^ +Po,e 


dsk 


qid = 0, s > 0, Vd, 


( 22 ) 


where [Q]d is the df^ column of Q and u \s a, row vector with = J2^=i +Po,£- This 

implies that lyQ = 0, and since Q is a generator we conclude that 

K 


T^d = 


E 

k=l 


gk- 


dRdjs) 

dsk 


+ P0,d — T^dX, 


where x does not depend on d. Observe now that by Lemma l5.11 we have 

i, ^ \ ~ dRd{s) 

Pd{s)-po,d = 2^gk 

k=i 
— TTf/X 

= H^{Z=d}]x - Po,d. 


Since Pd(s) = E 




, this implies that x = E 




= P{s). This 


shows that the environment becomes independent from the joint queue-length process in the 
heavy-traffic limit. □ 
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The flow equation, Eqn. (I14p . simplifies considerably in heavy traffic, as shown in the fol¬ 
lowing lemma. 

Lemma 5.3. If exists, then Rd{s) satisfies the following equation: 


0 = y^^Fk,d{s) 

k,d 


dRdjs) 

dsk 


Vs > 0, 


with Ffc^rf(s) defined as 


I^k,d{^) ■— 9k j ^ ^ hk,dSk 


Proof. We start by multiplying through Eqn. (121|) with N, followed by summing over d. Due 
to Q being a generator, this eliminates the right hand side with the transition rates qid- 


k,d 




E 




i=i 


dsi 


k^d _ 

Taking the limit iV —>■ oo yields 


hk,d9k 1 'Sfc T 


N J ds 


{N) 

d (g-s/Af) 


PoJ 
+ 0{N-^) 


= 0 . 


g. _ dR(i{s) dRd{s) 

0 ^ ^k,d^k / ^ 9j ^ fdk^ddk^k " 


k^d 


k,d 


dRd{^) sr^ \ sr^ 

/ ,9k o 2^'^j,dSj / y 9k,d9kSk 


dSh 

k,d j k,d 


dS}^ 

dRdjs) 

dsk 


^ ^ 9k j ^ ^ ^j,d^j h-kjd^k 
k,d 


dRdjs) 

dsk 


y~l Fk,d[t 

k,d 


dRdjs) 

dsk 


In what follows we focus on 


P'k,ooi^) • ^ ^ lP'k,d{s)TTd 9k 1 ^ ^ ^j,ooSj 9‘k,ooSk j i 


and denote its vector counterpart by Foojs) = {Fi^oojs),... ,Fk,oo{s)). 


(23) 

□ 
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6 Queue length distribution in heavy traffic 

We now state and consequently prove our main result about the queue length distribution. 


Theorem 6.1. When scaled by 1/N = (1 — the queue-length vector converges in 

distribution as (Pi^, • • • , P^k,L) iPi,oo, • • •, Pk,oo) i-e., ^ 1, namely 




(Mr’,..., Mf)) • l^z=d} A (Ml,..., Mx) • l^z=d} = ^d 


N 


Pl,oo pK,oo 

9l ' ' 9K 


X (24) 


where —>■ denotes convergence in distribution and X is exponentially distributed with mean 

Pk,oolpk Pd) 


EX = 


Coo Xyfc Pk,oo/idkPk) 
with pd := c'^^ Yhk ^k,d/pk and a = (oi,... , 00 )'^ being a solution of 

[Q ■ a]d = Cd{l - Pd)- 


(25) 


We will prove this theorem in the two following subsections, first showing in Section [6.II the 
state-space collapse observed in Eqn. (1241) and then in Section 16.21 we will show that X is 
exponentially distributed with the mean given by Eqn. (1251) . 


6.1 State-space collapse 

The first part of the proof of Theorem 16.11 is the state-space collapse. In this section, we 
assume limiv-j-oo lrf^^(e“®A) exists. 

Observe that due to Proposition 15.21 


Rdis) = E 


= E 


1 _ 

Y.k 9kMk 

1 _ Q-J2k^kMk 


HT,kMk> 0 } ■ Hz=d} 


■ 1 , 


EkdkMk i^kMk>0} 
=■■ R{s)7rd, 


■ TTd 


(26) 


where the last equation defines R{s), which is independent of d. We now derive some prop¬ 
erties of R{s). 

Lemma 6.2. R{s) is constant on a {K — 1)-dimensional hyperplane T-Lc, where 

"He := I s > 0 : ^ Sk = c| , c > 0. 


9k 
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Proof. We follow closely the steps of the proof of Lemma 3 in [33]. The proof has 3 steps: 
(i) Show that oo(s) is parallel to the hyperplane. Hence, any flow corresponding to oo 
that starts in the plane, stays in the plane, (ii) Show that R{s) is constant along each flow in 
the hyperplane and (hi) show that each flow in the hyperplane converges to a unique point. 
This implies that R{s) is constant on the hyperplane. 


(i) Foo{s) is parallel to Rc 

Observe that with 1 = '^kPk,oo and pk,oo = >^k,oo/Pk,oo, 




^ ^ Pk^oo I ^ ^ ^j,oo^j Pk,oo^k 




Pk,ooP'k,oo^k 


,oo^k E ,c 


^Sk 


= 0 . 


This indicates that the iL-dimensional vector Foo{s) is parallel to the hyperplane. 


(ii) R{s) is constant along flows in Re 

For each state s > 0, there exists a unique flow f{u) = (/i(u),..., fK{u))'^ parametrized by 


u > 0, such that 


/(O) = s and = Fk,oo{f{u)). 

Due to (i), any flow that starts in Rc, stays in Re- Now, 


(27) 


K 


dR{f{u)) dfk{u) dR{s) 


du 


du dsk 


E 

k=l 

^ ^ \\ dR(s) 

J^i^fc,oo(/(«))-^^ 

k=l ^ 


s=f{u) 


s=f{u) 


, .. dR(s) 

k=l d=l ^ 


s=f{u) 


K D 

k=l d=l 

= 0, by Eqn 


dRd(s) 


dsi 


s=f(u) 


implying that R{f{u)) is constant along each flow f{u) which lies in Rc- 
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(iii) Each flow in Tic converges to a unique point 

Here we first write the flow specifications in a vector-matrix form, then show that one eigen¬ 
value of that matrix is zero with eigenvector s* G "Hi, and the other eigenvalues are negative, 
and thus we can write f{u) = c • s* -|- g{u) where hmti_j.oo fl'(u) = 0. 

Eqn. (j27p can be written in matrix-vector form as 


f{u) = Af{u), 


with 


A = 


^ ffl(-^l,oo hl,oo) 51'^2 ,oo 

52^1,00 92 (^ 2,00 — h2,oo) 


gi^K,c 

92 ^K,c 


\ 9 kM ,00 9 k{^k ,00 l^K^oo) / 

Let D be the diagonal matrix with di = Pi^ooj9i on the diagonal. Then with 


^51(^1,00 — ^1,00) 52!^ A2, 


S := DAD-^ = 


91 


9'^(.^2,00 / 1 ' 2 ,Oo) 


9K 

9K 


Pl,o 

PK,c 

P 2.0 

PK,c 


-^K,c 

-^K,C 


^ PK,oo \ 


9 k(^k ,00 I^K^og) j 


is a generator corresponding to a hnite-state Markov chain. 

From the proof of Lemma 4 in [33], it is easily seen that the Markov chain corresponding to 
is irreducible (since we assume that all Afc,^ > 0 for all k and at least one d). Retracing 
the arguments stated there, for completeness, it follows that this Markov chain has a unique 
equilibrium distribution (column) vector, rj, such that = 0. In particular, 0 is an 

eigenvalue with multiplicity one and all other eigenvalues have a strictly negative real part, 
see ID. Since the eigenvalues of and A are the same, 0 is also an eigenvalue of A with 
corresponding right eigenvector s* = s* > 0,s* G Tii- The solution of the linear 

system f'{u) = Af{u),f{0) G Tic can now be written as the sum of the homogeneous and 
the particular solution, i.e. f{u) = c ■ s* + g{u), where lim^i-^oo ^((m) = 0. This implies that 
all the flows in "He converge to one common point c- s*. 

Combining (i), (ii) and (iii), we conclude that the function R{s) is constant on Re- D 


As a consequence of Lemma 16.21 the function R{s) depends on s only through the sum 
'^k=i(Pk,oo/9k)sk- Therefore, there exists a function i?* : M ^ M such that R{s) = 
R*{Y.k=iiPk,oo/9k)sk)- Then 



Pfc,oodR*(u) 

9k du 


^ — 1 iPk,oo / Qk'j^k 


18 













so we obtain 


D 




Af—>00 ' 

d=l 

D K 

d=l k=l 

^ ^ dR{s) 

-VTd 

d=l k=l 
K 

^ ^ Pk,a 
k=l 

dR*{v) 


dRdjs) 

dsk 


dsk 
dR*{v) 


dz; 


— Efc— 1 (pfe,oo/5A:)'^fc 


dz; 


P — Efc —1 (Pfc,oo/ Qk^^k 


(28) 


which only depends on v = J2k=iiPk,oo/9k)sk- Since we also have 


E[e' 


-Ef=lSfe^fe] = ]E e 32 ^^(p2.oo^^ Pi,00^1) 


-P£^Sk( 

SK \PK,cx> PI.00 J 


this together with Eqn. (i^ implies that = 0, for all j = 1,... ,K. 

Thus [pk/Pk,oo)Mk = {gj/pj^ca)Mj, for all k,j, almost surely. Combining this finding with 
that of Eqn. ([22]), we obtain Eqn. (i2l]l with X distributed as {gi/pi,oc)Mi. 

Remark 3 (Continuously modulated service requirements). In Section [3] we saw that the 
critical load is indeed reached when poo —>■ Ij since then po,d 0- This indicates that (1 —poo) 
is the right heavy-traffic scaling when pLk,d = Pk^d- This is less clear for a general Pk^d-, that 
is, for continuously modulated service requirements, where the environment can influence the 
departure rate of customers present in the system. For that setting, the workload process 
is no longer independent of the employed scheduling discipline, since the decision on which 
class to serve impacts the rate at which customers leave. We are not aware of any results on 
workload and waiting time distributions where the service distribution is a general function 
of both class and environment. 

The majority of the preceding queue length results in this paper can however be proven 
without the restriction of the product form, i.e., for continuously modulated service require¬ 
ments. The traffic intensity for this variant is defined as for the multi-class model above, 
only this time one cannot split the average class-A: service rate into pk,oo = PkCoo- The traffic 
intensity per class k, pk^oo = ^k,oo/Pk.oo-, is in line with the Markov-modulated single-server 
queues. Assuming there exists a scaling f{N) such that f{N){Mi ,..., MK)'i-{z=d} converges 
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in distribution, it can be shown that the empty probabilities po.d vanish in heavy traffic as 
N oo, for a general fik,d- This property then follows from Proposition 15.21 without relying 
on the workload results from Section [3] and the product form assumed there. Furthermore, 
under this assumption, all results in Section 16.11 hold, implying that a state-space collapse 
will appear. In other words, we can prove the first half of Theorem 16.11 However, we do not 
know what the distribution of the common factor X will be. 


6.2 Distribution of the common factor 


In order to prove that the limiting queue length distribution exists and to find the common 
factor of the queue length distribution in heavy traffic, the random variable X, we make use 
of the results on the workload of the total system. From [33] and Eqn. (I24D we have that 

w^y^ = x-y^. ( 29 ) 

In order to apply the workload result of Section [H] we first derive the service requirement 
of an arbitrary customer while being in state d. If Hk{-) is the distribution function of a 
class-A; customer’s service requirement, then the probability of a class-fe customer arriving 
and requiring service not exceeding x is ak,dHk{x). Summing over k now yields the desired 
distribution, 

K 

Hd{x) := '^ak,dHk{x). 
k=l 

The overall service requirement distribution thus depends on the state of the environment at 
its arrival. With exponential service requirements, the corresponding LST is given by 


K 


hd{s) = Yj 


^k,df^k 

and the first and second moment are given by 

ak,d 


, s > 0, 


hdi = Y—’ h,, = 2Y 

I I. 7. 




o^k4 


(30) 


(31) 


We can now apply the result of Theorem 13.31 with moments as given in Eqn. (|3ip . Hence, 
we have that W is exponentially distributed with mean 


EW = c: 


-1 


^ ^ Pk,oo/k-k ^ ^ C(j7r|jQ|j(l Prf) 
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where a is a solution of [Q ■ a]^ = — Xd Ylk "7^ ~ — Pd)- Along with Eqn. (f29]l this 

yields the mean of the exponential random variable X: 


EX 


ElE 

Xyfc Pk,oo/idkPk) 

Pk,oo/Pk X/d Pd) 

Coo Xyfc Pk,oo/{dkPk) 


(32) 


The hrst term of the numerator is in accordance with the results of |30] and [33], the second 
term is a result of the random environment. 

The results in Sections [51 and [6T] are based on the assumption that hm 7 v->-oo ' ^{z=d} 
exists. Since the scaled workload is tight, see Section jS] so is the scaled queue length. Then, 
by Prohorov’s theorem (|6|) there exists a subsequence of N such that ^Mk converges in dis¬ 
tribution, and hence for this subsequence limjv^-cxD exists. Since each converging 

subsequence yields the same limit, the limit itself exists (see corollary page 59 in [6]), i.e. 
i(M,Z = d) A M ■ l{z=d}-, as X —oo, with the limiting vector as in Eqn. (l2^ . 

This concludes the proof of Theorem 16.11 


7 Conclusion and future work 

We hrst studied the workload for a queue with modulated arrivals, service requirements 
and service capacity, and derived that the scaled workload converges to an exponentially 
distributed random variable in heavy traffic. The workload results obtained are valid for 
any service distribution and for any service discipline which does not depend on the envi¬ 
ronment. We then focussed on the special setting of a multi-class queue under the DPS 
policy and showed that the joint queue length distribution for such a system undergoes a 
state-space collapse in heavy traffic. Under the scaling of (1 —poo), the vector-valued limiting 
distribution is independent of the modulating environment and converges in distribution to 
a one-dimensional random variable times a deterministic vector. In this derivation, the dis¬ 
tribution of the scaled workload is a key quantity. With this we extend known results about 
the DPS queue to a Markov-modulated setting. 

Clearly an interesting question for future consideration is whether the state-space collapse for 
the DPS policy carries over to continuously modulated service requirements, as discussed in 
Remark 131 Another open question concerns the characterization of the moments of the queue 
lengths for the modulated DPS queue, outside of heavy traffic. Last but not least, modulating 
the weights of the DPS would open the possibility of dynamical scheduling based on the 
environment. The latter would be a study on its own, as already the stability conditions will 
no longer be independent of the weights of the DPS policy. 
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Appendix: Proof of Theorem 13.3 


The proof Theorem 13.31 is based on Theorem 4 in [T2] , which can be adapted to our model 
as follows: 

We start with notation and some preliminaries. Let A = diag(Ai,..., Xd), H{s) = diag(l — 
hi{s ),..., 1 — h£){s)), C = diag(ci,..., cd) and pg = (Po,i) • • • ■,Po,d)- Eurthermore Hi and 
H 2 are the diagonal matrices corresponding to the moments h^i and hd 2 , respectively, for 
d = 1,..., D. Recall Eqn. Q, [Q ■ a]^ = Cd — — Coo(l — Poo)- We will now construct a 

partial inverse of Q to make it easier to find a vector a which solves this equation. Let Qi 
and R be matrices such that 


Qi = 


<?22 d 23 


(l2D 


, R = 


\qD2 qns ■ ■ ■ qdd / 


/O 0 ••• o\ 

0 
: 

Vo / 


Then detQi 7 ^ 0 and due to Q being a generator (for more details see [l2]), we have 


QR = 


(Q ~^2 —TT3 . . . —TTp 

TTl TTl TTl 

01 0 ••• 0 

0 0 1 0 ••• 


VO 
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It follows that for any vector x, it holds that 

xQR = X — —TT. (33) 

TTi 

Then it can be verified with straightforward calculations that 

<2 — (oi, . . . , CLd') — R{ci Ai/iii Coo(1 Poo)) ■ ■ ■ ) CD XDhDl Cqo( 1 Poo)) 

= R[C - AHi]e - Coo(l - Poo)Re 

= R[C - AHi]e - r (34) 

is a possible solution vector, with r := Coo(l — Poo)R^- 

Define the vector (p = {ipi{s),... ,pd{s)) and write Eqn. ([5|) in matrix-vector terms, 

(p{s)Q = (f{s)[AH{s) - sC] + sp(^C. (35) 

Observe that, according to Eqn. © 


PqCs = Coo(l - Poo)- 


Now multiply from the right both sides of the new vector-matrix equation, Eqn. (1351) . with 
a D-dimensional vector of I’s, e, to obtain 

(p(s)[AD(s) - sC]e + sCoo(l - Poo) = 0. (36) 

Multiply from the right both sides of Eqn. (1351) with the matrix R to get 

ip{s)QR = (p(s) [AH (s) — sC]R + sp^CR. 

Rewrite this equation by using the property of Eqn. (l33]) to obtain 

<p{s) = -I- (p(s) [AH(s) — sC]R + spqCR. (37) 

VTi 

Iterate Eqn. (1371) with itself by inserting <p(s) into the right hand side of the equation to 
obtain, after some algebraic transformations, 

(^(s) = ^lM7r[/-FG(s)R]-hy(s) (38) 

TTi 


with G{s) := AH{s) — sC and 

y{s) := <p(s)[G(s)R]2 + sPqCR[G{s)R + /]. 
Substitute Eqn. (f38l) into Eqn. (l36)) to obtain an expression for <pi(s), 
>i(s) 


0 = 


TTi 

TTi 

TTi 


-7r[/ G(s)i?] -h y{s) 


G{s)e + sCoo(l - Poo) 


[7rG(s)e -h 'KG{s)RG{s)e\ -h y{s)G{s)e + scoo(l - Poo) 
[S2(s)+D3(s)]+Di(s) 


(39) 


(40) 
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with Bi{s) = y{s)G{s)e + sCoo(l — Pcxi), -^ 2 ( 5 ) = 7rG(s)e and -63(5) = 7 vG{s)RG{s)e. 


The next step is to insert the scaling s s/N for each term. Recall that using the heavy 
traffic parametrization introduced in Section [21 we have (1 — poo) = 1/A^- Now observe that, 
as N ^ 00, 


H{s/N) 

s/N 


Hi, 


His/N - H {s/N) 


Therefore the limit 


G{s/N) 

s/N 


{s/NY^ 


AHi - G, 


H 2 
2 ■ 


is a constant and since |(p(s/A^)| < 1 and ^ 0 (see Eqn. ([7])), we have 

+ pi^^CR[G{s/N)R + /] 


s/N 


s/N 


¥>{s/N) (^) 
^0, 


G{s/N) 

s/N 


R 


+ p^o^'^GR[G{s/N)R + I] 


as N ^ 00. Combining the above we obtain 

Bi{s/N) = y{s/N)G{s/N)e + sCoo(l - Poc)/N 

= sCoo(l - Poc)/N + o{N-^) = sCoo/N^ + o{N-^). 

Then 


B 2 {s/N) = TT [KH {s/N) - sC/N] e 

= ttA [H {s/N) - His/N] e + tt [AHis/N - sG/N] e 

,H{s/N)- His/N, ,^^,2 
= ttA- (^s/N)"^ - {s/Nfe - scoo(l - Poo)/N 

= -7TA^{s/N)^e + o{N~^) - sCoo/N'^ 

D 

= -(s/N)^ 7rdArf/irf2/2 - sCoo/N"^ + o{N-^) 
d=i 

since 7rA.ffie = Y,d=i T^d^^dhdi = PooCoo- Furthermore, 


Bz{s/N) = TT [AH {s/N) - Cs/N] R [AH {s/N) - Gs/N] e 


[a" 7'!'-cl 

R- 


[ {s/N) \ 


[ {s/N) \ 


= 7v{s/N)^ [AHi -G]R [AHi - C] e + o{N-^) 

= -7v{s/N)^ [AHi -G]{a + r) + o{N-^) 

= - ^ TTrf [{ad + rd){Xdhdi - Cd)] {s/N)^ + o{N~^) 

d 
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due to ii [C — AHi\ e = a + r, see Eqn. Under the heavy-traffic scaling, 

r = Coo(l - Poo)Re = CooN~^Re 

is an o(l) term. Observe that 

-{B2{s/N)+B^{s/N)) 

D 

= (s/Nf '^d[>^dhd 2 /‘^ + {ad + o{l)){Xdhdi - Cd)] -b sCoo/N'^ + o{N~'^). 


d=l 


Rearranging Eqn. (|4Up yields 


^Pl{s/N) = TTi 


Bi{s/N) 


= VTl- 


= TTl- 


-{B2{s/N)+B3{s/N)) 

CooS/N"^ + o{N~‘^) 


s/Ny '^d{Xdhd 2 /‘i‘ + adiXdhdi - Cd)] + Coos/N"^ -b o{N 2) 

1 + 0 ( 1 ) 


1 + Coo^ Ed T^d [Ad/id2/2 + ad{Xdhdi - Cd)] s + o(l) 

Let M be the desired mean stated in Theorem 3.3, that is 

A'l ■— Cqq ^ ^ T^d Xdhd2/‘^ + ad{Xdhdl Cd) 
d 

Then, taking the heavy-traffic limit, 

lim ifi{s/N) = lim ^ 

N—^oo N^oo TTi 1 + Ads 

i.e. the LST (f{s) converges in distribution to the LST of an exponentially distributed random 
variable with mean Ad. 
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